Numerical evidence against a conjecture on the cover time of planar graphs 
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We investigate a conjecture on the cover times of planar graphs by means of large Monte Carlo simulations. 
The conjecture states that the cover time t(Gjv) of a planar graph Gn of A'^ vertices and maximal degree d 
is lower bounded by t(Gjv) ^ CdN{\nN)^ with Cd = (ci/47r) tan(7r/d), with equality holding for some 
geometries. We tested this conjecture on the regular honeycomb (d — 3), regular square (d = 4), regular 
elongated triangular (d = 5), and regular triangular (d = 6) lattices, as well as on the nonregular Union Jack 
lattice (dmin ~ 4, dmax = 8). Indeed, the Monte Carlo data suggest that the rigorous lower bound may hold 
as an equality for most of these lattices, with an interesting issue in the case of the Union Jack lattice. The data 
for the honeycomb lattice, however, violates the bound with the conjectured constant. The empirical probability 
distribution function of the cover time for the square lattice is also briefly presented, since very little is known 
about cover time probability distribution functions in general. 
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I. INTRODUCTION 

The cover time of a graph is a classic problem in theoretical 
computer science and graph theory with many practical impli- 
cations, e.g., in the development of query processing and rout- 
ing algorithms in computer networks and distributed systems, 
and has attracted the attention of computer scientists, mathe- 
maticians, and physicists for more than 30 years [1-22]. 

For a finite, connected graph Gn = i^jE) of order N, 
the cover time t{Gn) is the maximum expected time over the 
possible starting vertices v E V it takes for a random walker 
jumping through the edges of Gn with uniform probabilities 
to visit every vertex of Gat at least once. Exact expressions 
for the cover time are rare except for the simpler graphs, e.g., 
for the complete graph, for which the problem reduces to the 
well-known coupon collector's problem, and for the path, cy- 
cle, and star graphs, among a few others [10, 15, 20, 21]. 

Most results on graph cover times come in the form of 
bounds on their expectation values, although little is known 
about the limit distributions of the related quantities. The exis- 
tence of a lower bound follows from t{Gn) ^ N — although 
it can be proved that ¥{t{Gn) ^ cN) ^ e""^, with a > 
depending only on c and (imax(GAr), the maximal degree of 
the graph [23] — , while the existence of an upper bound fol- 
lows from the recurrence of the associated Markov chains. 

For planar graphs, Jonasson and Schramm showed that [17] 

limiiif-JJ^^^Q, (1) 

where Gd > are universal constants depending only on 
d ^ d,„ax(GAr). The functional form T(GAr) - N{\nN)'^ 
for the cover time of the square lattice had been guessed ear- 
lier on the basis of Monte Carlo simulations and scaling anal- 
ysis, where a multiplicative correction (1 + c/\nN) to this 
form was detected, with c (the magnitude of the leading scal- 
ing correction) a constant depending on the boundary con- 
ditions of the finite graphs [5, 6, 9]. Following a sophisti- 
cated probabilistic-geometric analysis of the cover time of the 
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square lattice, Dembo et at. [19] conjectured that for the hon- 
eycomb, square, and triangular lattices, corresponding respec- 
tively to d — dmax — 3, 4, and 6, the constants Gd appearing 
in Eq. (1) are exactly given by 

C, = /tan(';), (2) 

further conjecturing that, with this Gd, inequality (1) may ac- 
tually hold as an equality for these geometries. 

In this article we investigate numerically the above men- 
tioned conjecture on some planar graphs to check whether the 
bound (1) with the conjectured constant Eq. (2) holds as a 
lower bound or as an equality. As we shall see, the empir- 
ical data for the cover time of the honeycomb lattice indi- 
cate a violation of the lower bound with the constant given 
by Eq. (2), requiring a smaller constant, while for some other 
planar graphs the conjecture seems to hold valid. 

II. MONTE CARLO DATA 

We investigate the conjecture encoded in Eqs. (l)-(2) on the 
regular honeycomb {d = 3), regular square {d = 4), regular 
elongated triangular (d = 5), and regular triangular {d = 6) 
lattices, and also on the nonregular Union Jack lattice ((i,„in = 
4, dmax = 8). These lattices are depicted in Figure 1. 

We computed the cover times on graphs with N — L x L 
vertices under periodic boundary conditions, with 256 ^ £ ^ 
1536, i.e., on graphs with up to 2.359 x 10^ vertices. For each 
graph geometry and size, t{Gn) is obtained as an average 
over 10^ samples. Our Monte Carlo data together with the 
conjectured values t*{Gn) = GdN{lnN)'^ obtained from 
Eqs. (l)-(2) appear in Figure 2. 

Notice that the toroidal graphs obtained under periodic 
boundary conditions are not planar, although they are locally 
very close to planar. The really important fact in a planar 
graph for the cover time problem, however, is that its edges do 
not cross, not that it can be embedded in a plane. Moreover, 
finite graphs with open boundary conditions cannot be regu- 
lar, since the vertices at the boundaries are of a smaller degree. 
The asymptotics in the two cases (open and periodic boundary 
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conditions) are expected to be the same, and most results on 
graph cover times, including those to which we want to com- 
pare our own results, are obtained for graphs under periodic 
boundary conditions. 

According to the conjecture, we must observe 
T* {Gn)/t{Gn) ^ 1 for all planar graphs, a condition 
that our numerical data does not support for the honeycomb 
lattice. Other lattices observe the bound, with the square 
lattice being "borderline." A naive extrapolation of the ratios 
T* {Gn)/t{Gn) would give an extrapolated value greater 
than 1 in almost all cases depicted in Fig. 2. This indicates 
that possible corrections to t*{Gn) must go in the direction 
of decreasing its value. Since previous results in the literature 
suggest just the opposite, i.e., that, if anything, t*{Gn) may 
be missing a (1 + c/\nN) correction with c > [5, 6, 9], 
we are led to believe that the constant with d — (i,„ax 
is overshooting. Supplementary evidence comes from the 
behaviour of the Union Jack lattice with respect to d. While 
the bound (1) requires d ^ rfmax(GAr), our data suggest 
that this requirement is probably not optimal. We plot the 
ratios t*{Gn)/t{Gn) for the Union Jack lattice both with 

d = dmax = 8 and with d = d = ^(dinin + ^max) = 6, 

the average degree of the lattice, and we found that the ratio 
with d ~ d provides a better lower bound than the ratio with 
d = dmax; see Fig. 2. This makes us wonder if the average 
vertex degree 

' ' vev(G) 

where d{v) is the degree of vertex ti, is not a better constant to 
be used on purportedly uruversal formulas for planar graphs 
than the maximum degree dmax- This could be tested on pla- 
nar random graphs — e.g., on Delaunay triangulations of ran- 
dom points on the plane [24] — , for which d can assume dif- 
ferent values, integer or not. Recently, planar random graphs 
have attracted the attention of physicists and mathematicians 
interested in their connectivity and percolation properties [25- 
30], but their cover times remain unexamined. 

III. THE EMPIRICAL PROBABILITY DISTRIBUTION 
FUNCTION OF THE COVER TIMES OF THE SQUARE 
LATTICE 

Very little is known about the probability distribution func- 
tion (PDF) of t{Gn)- Actually, it is an open problem to prove 



mean 


m 


107111924 


variance 


s 


12 308 707 


skewness 


91 


0.91374 


excess kurtosis 


92 


1.51095 



TABLE I. First moments of the empirical probability distribution 
function of the cover time of the square lattice depicted in Fig. 3. 



that t{Gn) has a nondegenerate limit law [15, 19, 20]. It 
seems that the only result on this regard available to date is 
a concentration result that states that, under mild conditions, 
the cover time is well approximated by its expected value as 
N oo [1]. It is thus of some interest to explore our empir- 
ical data to characterize the PDF of the cover time, although 
we will not attempt to identify or infer this distribution here. 

Figure 3 shows the histogram plot of 10^ cover times sam- 
pled for a square lattice of = 1280 x 1280 vertices. For 
these data, we compute the sample mean rn = (Ti(G'Ar)) and 

first few central moments mk = {{Ti{GN) — 'm)^)^ from 
which we compute the sample standard deviation s = ^m^, 

skewness (7i = tos/toj^^' ™d excess kurtosis 52 = m^jrn^ — 
3 [31]. The values of these quantities are collected in Table I. 
We also fitted the data to the beta PDF given by 

rescaled to the interval [a, 6] with a = min{Ti(G'Ar)} ~ 
73 391036 and 6 = max{n(Gjv)} = 213 039197. We 
choose the beta distribution because it has finite support and 
can take several different shapes; moreover, beta-like PDFs 
for the cover times of some special graphs were found in pre- 
vious investigations [5, 6, 21]. The empirical data together 
with the adjusted beta PDF appear in Fig. 3 [32]. 

The positive skewness g\ indicates that the empirical PDF 
is right-tailed, with the bulk of the observed values lying to the 
left of the mean, although this feature is not very clear from 
the histogram (3) because gi is not very large. The moderately 
high value of the excess kurtosis 52, in turn, indicates that the 
empirical PDF is markedly non-normal, with a sizeable pro- 
portion of the data in its right tail contributing to the variance 
observed. We notice that after 10® samplings, the ratio sjm 
has stationed at ^ 11.5%. 
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no. 2. (Color online) Cover times for the planar lattices depicted in 
Fig. 1. Each graph displays the empirical cover times t(Gjv) (open 
squares, left scale) together with the conjectured value r*(Gjv) = 
CdN{lnN)^ (solid line, left scale) and the ratio t*{Gn)/t{Gn) 
(full squares, right scale). Each point of the t(Gn) curves was ob- 
tained as an average over lO'' samples. For the Union Jack lattice we 
display the ratio r* (Gjv)/r(Gjv) both for Cd with d — dmax = 8 
and d = d = |(dmiii + rfmax) = 6; the plots of t*(Gn) for these 
two values of d are indistinguishable on the left scale of the graph 
and appear as a single line passing through the open squares. 




FIG. 3. (Color online) Histogram plot of 10 cover times sampled 
for a square lattice of N — 1280 x 1280 vertices. There are 23 
bins in the histogram, with the leftmost bin centered at 70 x 10® and 
the righmost bin centered at 180 x 10®. The vertical dashed lines 
indicate the empirical mean m and ±s intervals. The continuous 
line corresponds to the adjusted beta PDF (4) with shape parameters 
a = 5.759±0.008and/? = 18.04±0.02. The empirical cumulative 
distribution function is also shown (open squares, right scale). 



IV. CONCLUSIONS 

Our numerical data for the cover time of the honeycomb lat- 
tice provide evidence against the conjecture set forth in [19] 
regarding the constant Cd appearing on the lower bound (1) 
for planar graphs. Otherwise, for the other lattices investi- 
gated in this article the functional form given by Eqs. (l)-(2) 
seems to hold valid, possibly as an equality. In summary, for 
d = 3 our data seem to falsify the conjecture, for d = 4 it is 
"borderline," and for d ^ 5 the conjecture holds easily. That 
there must be something special about the d = 4 case has 
been long recognized [15], and this might have showed in our 
finite-size simulations. Notice that the lower bound (I) re- 
mains valid for some Cd, just not with the Cd given by Eq. (2) 
with d = dinax- Since the Cd given by Eq. (2) is monotone 
decresing in d, it can be used in Eq. ( 1 ) to validate the bound, 
but then necessarily with some d > d,nax- 

It can be argued that our data were taken too short from the 
limit N oo. However, if the conjecture is to be saved, 
that would mean a nonmonotone convergence of the ratios 
r* {Gn)/t{Gn), something that our data do not indicate. 

Any case for the average degree d given by Eq. (3) would be 
welcome. It may be that the situation here is similar to that of 
the determination of the critical percolation threshold pc, for 
which good approximations and scaling relations were found 
based on the mean Euler characteristic of the critical percola- 
tion patterns [29]. It should be remarked, however, that for the 
critical percolation threshold problem one does not expect to 
find a "universal" formula in terms of the maximal vertex de- 
gree d alone, simply because of the empirical observation that 
several different lattices with the same d have disparate pc- 
Notice that the constant Cd in Eq. (2) is closely related with 
the filling factor / introduced by Suding and Ziff in order to 
relate to the number of sites per unit area of Archimedean 
lattices [33], 
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where (n^^rij^, . . .) is the Griinbaum-Shephard representa- 
tion of the Archimedean lattice — e.g., the square lattice is de- 
noted by (4^), while the elongated triangular lattice is denoted 
by (3"^, 4^) [34]. The investigation of different planar graphs 
of same degree (regular or average) — e.g., for d — 4, the 
square, kagome, and Archimedean (3,4,6,4) lattices — may 
help to elucidate the above questions of geometric character. 

Finally, our exploratory analysis of the PDF of the cover 
times of the square lattice in Sec. Ill is admittedly jejune; that 
was not the focus of this work. We were nevertheless able 
to establish that the empirically observed PDF is not normal, 
with a leptokurtic shape (52 > 0) . A proper investigation of 
the PDF of the cover times of planar graphs — e.g., by model 



selection among candidate left and right limited, univariate 
distributions — is still lacking and provides an interesting di- 
rection for further statistical work on cover times. 
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